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We investigate consensus formation and the asymptotic consensus times in stylized individual- 
' or agent-based models, in which global agreement is achieved through pairwise negotiations with 

or without a bias. Considering a class of individual-based models on finite complete graphs, we 
5_l . introduce a coarse-graining approach (lumping microscopic variables into macrostates) to analyze 

the ordering dynamics in an associated random-walk framework. Within this framework, yielding a 
linear system, we derive general equations for the expected consensus time and the expected time 
spent in each macro-state. Further, we present the asymptotic solutions of the 2- word naming game, 
and separately discuss its behavior under the influence of an external field and with the introduction 
C^l ' of committed agents. 

r Z?- PACS numbers: 87.23.Ge 05.40.Fb 

q , Individual- or agent-based models provide invaluable tools to investigate the collective behavior and 
Q ■ response of complex social systems. These systems typically consists of a large number of individuals 
interacting through a random and sparse network topology. Despite this non-trivial network topology, 
it has been demonstrated in many recent examples that unlike in low-dimensional spatially-embedded 
systems with short-range connections, the collective dynamics in sparse random graphs (with no 
community structure) exhibit scaling properties very similar to those observed on the complete graph. 
Therefore, studying fundamental agreement processes on the complete graph can yield insights for 
the ordering process in more realistic sparse random networks. In this paper we consider two simple 
individual- based models and develop a mathematical framework which yields asymptotically exact 
consensus times for large but finite complete graphs of size N. In particular, after demonstrating the 
feasibility of this framework on known examples, we apply it to study two distinct stylized approaches in 
social influencing: (i) influencing individuals by a global external field (mimicking mass media effects) 
. and (ii) introducing committed individuals with a fixed designated opinion (who can influence others 
\^ ' but themselves are immune to influence). In the former case, we find the external field dominates the 
T^j- , consensus in the large network-size limit, while in the latter case we find the existence of a tipping 
■ point, associated with the disappearance of the metastable state in the opinion space. The results 
further our understanding of timescales associated with reaching consensus in social networks. 
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I. INTRODUCTION 

Research on non-equilibrium models for social and opinion dynamics has attracted considerable attention [1-3]. 
In this paper we present a method to calculate consensus (or ordering) times for a large class of individual-based 
models with discrete state-variables on the complete graph of size N. The mode of reaching consensus depends on 
the initial configuration and whether or not the individuals are under the influence of a bias. For example, starting 
from a configuration of evenly mixed opinions with no bias, these systems can exhibit "coarsening" in which case the 
consensus time typically diverges with the system size in a power-law fashion [3-9] . On the other hand, initializing the 
system in one particular opinion and exerting some weak bias favoring another opinion create systems that undergo 
an "escape" from the meta-stable state [10]. For networks of finite size, it is important to characterize the consensus 
time, i.e. the time until the system fully orders. 

Direct simulations of the above behaviors can be time-consuming for large systems and essentially impossible even 
for moderately-sized systems when the system initially is in a meta-stable configuration, as the escape time can increase 
exponentially with the system size. The method presented here provides a way to obtain the asymptotic behavior of 
consensus times, including the cases associated with extremely slow meta-stable escapes. While the method and the 
results presented here are applicable to the complete graph, consensus times often exhibit the same asymptotic scaling 
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TABLE I: Update events for the voter model and the associated random walk transition probabilities 
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with A?" in large homogeneous sparse random networks [9, 11, 12], hence one can gain some insight how ordering and 
consensus evolves in realistic social networks. 

The Naming Game (NG) first emerged in linguistic modeling of the transition from microscopic local interactions 
to global consensus in the absence of a global coordinator [3, 8, 13]. This phenomena of total synchrony or global 
consensus has general interest and applications such as sensor networks, some aspects of neuro-science models for 
brain functions and social dynamics. Earlier research in this domain mainly discussed two aspects of this model. One 
focuses on the relationship between the network topology and the dynamics [12, 14, 15], the other, in contrast, ignores 
all the effects of network topology and studies the intrinsic property of the dynamics on the complete graph [8, 16]. 
The latter for large network size N (also referred to as mean-field or homogeneous mixing) is what we consider in this 
paper. We assume here an initial configuration where every agent has a non-empty name list of words. The dynamics 
on the complete graph is driven by the mean field (in this context referring to the fraction of nodes in each state) 
rather than the detailed network states (including the states of all the nodes). It naturally leads to a coarse-graining 
approach by lumping microscopic variables into macro-states. Our paper studies the impact of a global external field 
(mimicking mass media effects) and of presence of committed agents (who will firmly stick to and convey a designated 
opinion) on the consensus process. 

A point worth mentioning is that the Naming Game constitutes an opinion dynamics model in which in which a 
node can possess multiple opinions simultaneously. Such models with intermediate states have only recently begun 
receiving attention [17, 18] in statistical physics literature, and we believe that this study is an important contribution 
to this body of literature. 

The paper is organized as follows. In Sees. II and III, we introduce the coarse-graining approach to map models 
for opinion dynamics on the complete graph to an associated random walk problem, and test our framework through 
known results for the spontaneous agreement process (without any influencing) [11, 16, 19-21]. We derive the equations 
in a linear-system form for the expected total consensus time and the expected time spent in each macro-state which 
can be solved in a closed form for the voter model [1, 4, 22] (Sec. II). Then, in Sec. Ill, we employ the same method for 
a variation of the Naming Game (NG) with two words [21]. Finally, and most importantly, we investigate models for 
social influencing: we present new asymptotic results and discuss the behavior of the 2-word NG under the influence 
of an external field (Sec. IV) and in the presence of committed agents (Sec. V). A brief summary is given in Sec. VI. 



II. CONSENSUS TIME IN THE VOTER MODEL ON THE COMPLETE GRAPH 

First, we consider a well studied prototypical model for opinion formation, the voter model [1, 4, 22]. In this model, 
the evolution of a suitably defined global variable can be easily mapped onto a random-walk problem. Further, the 
solution of this model is known in all dimensions, including the complete graph, hence our method can be tested. 

Given a network of N nodes, with each node in a state chosen from the set of possible opinions X, the voter model 
is defined by the following update rule: 

A pair of nodes, a "speaker" and a "listener", are chosen at random. The listener then changes its state to that of 
the speaker. 

If the set of nodes in the network is denoted by S, {\S\), then the above rule defines a Markov chain in an N 
dimensional space X s . Under mean field assumption, which is justifiable when dealing with complete graphs, one 
coarse-graining approach is to take all network states in X s corresponding to the same mean field n — (ni, n,2, Um) 
as a macrostate, where rij denotes the number of nodes in state i £ X and M = \S\. Therefore, the coarse-grained 
random process is valued in a M — 1 hyper-plane (since n i = N) in M dimension space. When M = 2 (2-state 

voter model), taking X = {A, B}, the coarse-grained process is on a segment (since ua + ub = N(nA,nB >= 0)), 
and all macrostates can be represented by a single discrete variable ua = 0, 1, N. 

In each time step An a, the change of ha, is a random variable depends only on current macrostate ha- Its possible 
values and the corresponding events and probabilities are listed Table I. 
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Hence, the coarse-grained process of the voter model can be mapped to a random walk in 1-d, 

T-l 

n A {T) = n A (0) + &n A (n A (t)) , 



(1) 



t=o 



where n A — 0,7V are two absorbing states. 



A. First-step analysis of consensus times in the voter model 



The expected time before absorption which is the quantity of the most interest can be evaluated by first-step 
analysis [23]. The idea is based on a straight forward statement: The absorption time can be decomposed into two 
parts, the time steps before and after leaving the current macrostate. The former one is called the residence time 
t r (n A ) of the given macrostate n A . We denote the expected time before absorption and the expected residence time 
starting from a specific macrostate n A on complete graphs with TV nodes as r(n A \N) and t(n A \N), respectively, or 
simply r(n A ) and t(n A ) when it is not ambiguous. Then taking expectation of all random variables mentioned in the 
statement above, we get the following equations: 



'(n A ) 



_E[time before leaving n A ] + E[time before absorption and after leaving n A ] 

E[t r {n A )\ + E[E[t\me before absorption and after leaving reentering new macrostate n' A ]] 



t(n A ) + 



Ei=-i i P(An A (n A ) = i)r(n A + i)) 



t{n A ) 



l-P(An A {n A ) =0) 
where t(n A ) is given by following argument: 

t(n A ) = E^P^tia) - k)k = ^ =1 P(t r {n A ) > k) 



-(T(n A + 1) + t(ji a - I)) 



(2) 



Y,f =1 P{An A {n A ) = 0) 



fe-i 



1 



1 - P(An A (n A ) = 0) 2^(1-^)' 



(3) 



Defining f = (r(l), t(N - 1)) T , f = (i(l), t(N - 1)) T and using boundary conditions r(0) = r(7V) = 0, we 
rewrite Eq. (2) as a linear system: 



T = Qf+t 



(4) 



where 



Q 



l 

2 



(5) 



(JV-l)x(JV-l) 



r(n A \N) can be solved exactly from this equation in terms of n A and TV. Define 8 l N as a N-entries column vector in 
which only the ith entry is non-zero and has value 1. We have: 



r(i\N) = 



(s^fii-Qr't = u-t 



(6) 



where u = (ui,u 2 , uj, ...,un-i) is the solution of equation (/ — Q) 1 u = (S l N _ 1 ). Intuitively, uj can be understood 
as the average number of visits (assuming entering and leaving a macrostate as one visit) of macrostate n A = j before 
absorption and starting from macrostate n A — i. Moreover Tj = Ujtj is actually the average number of time steps 
(including self-repeating steps) spent on the macrostate n A — j before absorption. It is quite easy to show that 
r(n|7V) is monotonic with respect to N. Therefore, in order to obtain the large N behavior of t(ii a \N) we focus on 
the case when N is even and always assume n A (0) = N/2. In such cases, u = (1, 2, 3, N/2, ...3, 2, 1) T and 



N(N - DEK-'jf^, + N a N(N-l) l~ f + N 

ln(2)/V(/V — 1) + N . (7) 
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As is the convention for agent-based models in statistical physics, unit time is assumed to consist of N update events. 
Thus, following this convention, the normalized consensus time is t(N/2\N)/N is ln(2)JV to the leading order. This 
agrees with the scaling behavior obtained previously for the characteristic (relaxation) time of the voter model through 
a Fokker-Planck [19, 20] approach and through simulations [11]. . 

Figure 1 shows the comparison of the consensus time from numerical simulations (averaged over 50 runs for each 
iV) with the analytical results. 



700 



600 



500 



400 



300 



200 



100 







' 1 1 1 , 

/ 

* / ' 




* numerical simulation 
analytical result: log2*N 


/ * 












* / * 
* / * 




* /* 
* / 




*/ * 








rill 



100 



200 



300 



400 



500 
N 



600 



700 



800 



900 



1000 



FIG. 1: Consensus time for voter model on complete graph. The vertical axis represent the consensus time (normalized by N). 
The horizontal axis is the number of nodes in complete graph. Each star point is an average of 10 runs of numerical simulations 
of voter model and the solid straight line consists of the solutions of the linear equation for each N value. 



III. THE 2- WORD NAMING GAME 



The Naming Game (NG) [15, 16, 21] is somewhat more complicated than the voter model because for the given set 
of all possible opinions (words in the original NG) X (\X\ = M), the state of each node is a member of the power set 
2 X - the set of all subsets of X - rather than X itself. Moreover, the update rule is replaced by: 

A pair of nodes, a "speaker" and a "listener", are chosen at random. The speaker then randomly selects one 
word from her list and communicates it to the listener. If the listener already has this word (termed "successful 
communication"), she deletes all other words in her list (i.e., collapses her list to this most recently communicated 
word); If the listener does not have the word communicated by the speaker, she adds it to her list (hence, individuals 
can carry more than one word at a time). 

The slight difference between the update rule defined above and that of the original NG [15] is that here, upon 
"successful" communication, only the listener changes its state. This restriction eliminates steps of size 2 in the 
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TABLE II: Update events for the 2-word naming game and the associated random walk transition probabilities 



speaker 


listener 


event 


An(n A ,n B ) 


probabilty 


B or AB 


A 


A — ► AB 


(-1,0) 


P(A-) = n A (N -UA + n B )/2N 2 


A or AB 


AB 


AB ^ A 


(1,0) 


P(A+) = (N-n A - n B )(N + n A - n B )/2N 2 


A or AB 


B 


B ^ AB 


(o,-i) 


P(B-) = n B (N + n A - n B )/2N 2 


B or AB 


AB 


AB ^ B 


(0,1) 


P(B+) = {N -n A - n B )(N - n A + n B )/2N 2 


A, B or AB 


A or B 


unchanged 


(0,0) 


P = (n A + n B )/2N + (n A - n B ) 2 /2N 2 



associated RW model, making it easier to apply the method developed in Sec. II. A. Furthermore, we will consider the 
version of the NG with only two words [21]. 

The coarse-graining approach mentioned in Section II merges all network states labeled by the same vector n = 
(ni) 2 M_ 1 ,i 6 2 X \{0} into a macrostate. Here rn is the number of nodes in state i. The coarse-grained random process 
takes values in the 2 M — 2 hyper-plane: J2ie2 x \{i>} Ui ~ tnus we can ma P l ne coarse-grained process into a 2 M — 2 
dimension random walk. In the case of 2-word NG M = 2, so assuming X = {A, B}, n = [ji{A}i n {B}i n {A,B}) or 
(riA,nB,nAB) for short, where ha, tib, and uab are the number of individuals carrying word A, word B, or both, 
respectively. Since uab = N — ua — ns, we dump one redundant dimension and take the 2-d vector n = (uajTIb) 
to represent the macrostate. In each time step, the change of macrostate An have five possible values. For all these 
possible values of An, the corresponding events and probabilities are listed in Table II. 

In Table II, the event A — > AB, for example, means the listener node changes its state from A to AB. Analogous to 
the procedure followed for the voter model, we map the coarse-grained 2-word NG to a 2-d random walk, 

T-l 

n(T) = n(0) + An(n(t)) . (8) 
t=o 

Here the domain D of n is the set containing all integer grid points bounded by the lines ua — 0, 115 = and 
riA + riB = N while n = (0,N),(N,0) are two absorbing states. The expected drifts, E[An(n)] = (P(A+) — 
P(A-), P(B+) — P(B—)) are plotted in Fig. 2 for N = 10. As shown in Fig. 2, on average, n(t) will quickly go 
to a stable trajectory and then slowly converge to one of the consensus states. In other words, in contrast to an 
unbiased random walk as in the voter model, the 2- word NG is " attracted" to the consensus state after a spontaneous 
symmetry breaking. So it is reasonable to expect that the 2- word NG will achieve consensus much faster than the 
voter model, both starting from the ua = tlb = N/2 initial configuration on the complete graph. 



A. First-step analysis of consensus times in the 2-word Naming Game 

We now repeat the first-step analysis (developed in Sec. II), noting that the method is essentially independent of 
the number of dimensions. Assume r(n\N) and t((n)\N) are expected numbers of time steps before absorption and 
before leaving the current state starting at the macrostate n = (n^,ns), then we have: 

,^ ^ , E Je {(i.o),(-i.o),(o.i),(o.-i)} P ( Ar? (") =i) r(n-H) 

TW - t{n) + 1-P(A*(*) = (0,0)) 

P(A+)T(nA + l,n B )+P(A-)T(nA-l,n B ) + P(B+)T(nA,nB + l) + P(B-)T(nA,nB-l) , n , 

= *(") + T - p(o) (9) 

and 

1 1 

t[n) ~ l-P(An(n) = (0,0)) ~ 1 - (n A + n B )/2N - (n A - n B ) 2 /2iV 2 (10) 

Ordering all the macrostates n (except the two absorbing states) in a vector, and arranging T(ri),t(n) in the same 
order to get f and t whose dimensions are (N + 2)(N + l)/2 — 2, we write the Eq. (9) in the same linear system form 
t = Qf+t. Furthermore, taking 6% as a column vector where the only nonzero entry is 1 corresponding to n, we solve 
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FIG. 2: Vector field E[An(n)} of the random walk coarse-grained from 2-word naming game. Each vector is the expected drift 
of the random walk at macrostates n. The network is 10 nodes complete graph and the domain of random walk is the lower 
triangle of the square lattice 



expected number of visits to each macrostate through the equation (I — Q) T u = S^. The expected number of time 
steps spent on each macrostate T(riA,nB) is obtained by multiplying corresponding elements of u and t. Although 
in this case the matrix Q is not easy to write generally, it has a good property that sum of each row is < 1 and for 
some rows (those corresponding to n = (N — 1,0) and (0, N — 1)) the inequality is strict. Consequently the moduli 
of all the eigenvalues of Q is strictly less than 1 and (/ — Q) is invertible, therefore the existence and uniqueness of 
the solutions are guaranteed. 

We study the case of even ./V and unbiased initial state n(0) = (N/2, N/2), and compare the normalized consensus 
time t(N/2, N/2\N) obtained from numerical simulations against the solution obtained from the linear system. Figure 
3 shows that the normalized consensus time has an order 0(ln N) which is smaller than the order O(N) for the voter 
model. Note that the O(lniV) consensus time of the original 2-word NG with the above initial configuration has been 
previously found using simulations and a rate equation approach [16, 21]. 



Figure 4, depicting the expected number of time steps spent on each macrostate T(n J 4,ns) when N = 100 and 
n(0) = (50, 50) shows that the time before absorption is mainly spent on the macrostates near the consensus states. 



IV. THE 2- WORD NAMING GAME WITH EXTERNAL INFLUENCE 



The previous section focused on the time taken by the system to spontaneously reach consensus. A natural question 
that arises is how consensus can be sped up through an external influencing force such as mass media [24, 25], In this 
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FIG. 3: Consensus time (normalized by N) as a function of network size N for 2-word naming game on complete graph. Each 
star point is an average of 10 runs of numerical simulations of 2-word naming game and the solid straight line consists of the 
solutions of the linear equation for each TV. 



TABLE III: Update events for the 2-word naming game with central influence and the associated random walk transition 
probabilities 



speaker 


listener 


event 


An(n A ,n B ) 


probabilty 


B or AB 


A 


A -> AB 


(-1,0) 


P(A-) = n A {N -n A + n B )/2N 2 


A, AB or f 


AB 


AB ^ A 


(1,0) 


P(A+) = (1 - f)(N -riA- n B )(N + n A - n B )/2N 2 + f(N -n A - n B )/N 


A or AB 


B 


B -> AB 


(o,-i) 


P(B-) — n B (N + n A — n B )/2N 2 


B or AB 


AB 


AB -> B 


(0,1) 


P(B+) = (1 - /)(JV -n A - n B )(N - n A + n B )/2N 2 


A, B or AB 


A or B 


unchanged 


(0,0) 


Po = {n A + n B )/2N + {n A - n B ) 2 /2N 2 



section we study this case: the 2-word NG subject to an external field of magnitude / for which the update rule is 
defined as follows: in each time step, if the listener is in a mixed state i.e. AB, it will with probability / spontaneously 
change into state A and with probability 1 — f follow the original NG rule (Sec. III). All the differences between this 
case and the spontaneous case lie in the An(n) which are listed in Table III. 

In Fig. 5 we show how the vector field _E[An] which is intuitively the "drift" part of the course-graining random 
walk changes at different influence levels /'s. 



Following the first-step analysis, we can solve for the expected consensus time r and the expected number of time 
steps spent at each macrostate T(n J 4,ns) starting from any given macrostate. A specific solution of T(n^,ns) on a 
complete graph with N — 100 starting from macrostate (50, 50) is shown in Fig. 6. As shown, there are two peaks 
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FIG. 4: The expected time spent on each macrostate before consensus in the 2-word NG on a complete graph with TV = 100 
nodes. The vertical axis T[riA, tib) is the expected time that the random walk spends in macrostate (nA,ns) before consensus, 
starting from the (t),a(0), n_e(0)) = (50,50) initial macrostate. 



around the two consensus states, just as in the spontaneous case, although the peak near the consensus state which 
the external influence prefers (all A state) is much higher than the other one, even at a low influence level / = 0.05. 



A. First-step analysis of probability of consensus 

To better understand the effect of external influence, we apply the first-step analysis on the probability of reaching a 
specific consensus state. Defining Pa(^) as the probability of going to an all- A consensus starting from the macrostate 
n = (nA,n,B), PA(n) follows: 

P A (n(t)) = E[P A (n(t + 1))] = J2 P(An(n) = i) P A (n + i) (11) 

je{(i,o),(-i,o),(o,i),(o,-i),(o,o)} 

and satisfies the boundary conditions Pa{N, 0) = 1 and Pa(0,N) = 0. Ordering Pa(") of all macrostates (including 
the two consensus states) in a vector Pa, we rewrite the equations as Pa = QoPa- Qo is a square matrix of order 
(N + 2)(N + l)/2 and all its elements are given in Table III. 

In Fig. 7, we consider the NG on 100-node complete graph starting from the macrostate (no,N — no). For uq = 50 
and central influence / — (the left end of the black curve), it has equal probability of going to all A and all B 
consensus. When no < 50, it is more probable to go to all B consensus without central influence. However, with a 
biased central influence /, one can always convert the preference of the process to the opposite side - all A consensus. 
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FIG. 5: Several vector fields are drawn in this figure with different colors, each color show the drift of coarse-grained random 
walk E[An(n)] on a complete graph with N = 20 nodes at different central influence levels /. The length of each vector has 
been rescaled to its square root to avoid cluttering the whole graph. 

Furthermore, we show in Fig. 8 that the external influence becomes more powerful in forcing the network to a desired 
consensus state when network size N grows larger. Starting from an unbiased macrostate (N/2, N/2), the probability 
of going to all B consensus 1 — Pa(N/2, N/2) (which is against of the central influence) decays exponentially along 
with the network size N. So it is reasonable to expect that in real social network for which network size N is very 
large, a very slight biased central influence can strongly affect the social consensus. 



A second method of speeding up consensus is by introducing an intrinsic bias in the system, through a set of 
inflexible agents [10, 26] promoting a designated opinion. We refer to such individuals as committed agents [10]. 
Intuitively, the introduction of committed agents will break the symmetry of the original NG and will facilitate global 
consensus to the state adopted by the committed agents. In this section we provide asymptotic solutions of 2-word 
NG with committed agents 

Suppose that the number of committed agents is n q and all committed agents are in state A. The corresponding 
events in the NG with committed agents and the associated RW transition probabilities are summarized in Table IV. 
The equations f = Qf + t and (I — Q) T u = <5^ are derived in exactly the same way as in the non-committed case. 

From the dynamics of infinite systems with homogeneous mixing, one can expect [26, 27], that there exists a critical 
value q c of the committed fraction q — n q /N , above which consensus times drop dramatically. More specifically, for 
q < q c the phase space exhibits three fixed points: a "meta-stable" one, dominated by individuals in state B, an 
stable absorbing fixed point with all individuals in state A, and a "saddle" point separating them. Initializing the 
system in a configuration corresponding to macrostate (n^(0), ng(0)) = (n q ,N — n q ) (a small number of committed 
agents embedded among Bs), the system quickly relaxes to the meta-stable fixed point and stays there for times 
exponentially large with the system size (i.e., forever in an infinite system). As q^q c , the meta-stable fixed point and 



V. 



2- WORD NAMING GAME WITH COMMITTED AGENTS 
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FIG. 6: Expected time spent on each macrostate before consensus T(n A , Ub) in the 2- word NG on complete graph with N = 100 
with central influence / = 0.05, starting from a unbiased initial macrostate (n A (0), n B (0)) = (50,50). 

TABLE IV: Update events for the 2-word naming game with committed agents and the associated random walk transition 
probabilities 



speaker 


listener 


event 


An(n.4,rifl) 


probabilty 


B or AB 


A 


A ->• AB 


(-1,0) 


P(A-) = (n A - n q )(N - n A + n B )/2N 2 


A or AB 


AB 


AB^ A 


(1,0) 


P(A+) = (N -n A - n B )(N + n A - n B )/2N 2 


A or AB 


B 


B ^ AB 


(0,-1) 


P(B-) = n B (N + n A - n B )/2N 2 


B or AB 


AB 


AB -> B 


(0,1) 


P(B+) = (N-n A - n B ){N -n A + n B )/2N 2 


A, B or AB 


A or B 


unchanged 


(0,0) 


Po = (n A + n B )/2N + (n A - n B f /2N 2 +n q {N -n A + n B )/2N 2 



the saddle point merge, and become an unstable fixed point. For q > q c , regardless of the initial configuration, the 
system quickly relaxes to the all-A absorbing fixed point. Figure 9 confirms the above scenario. Figure 9(a) and (b) 
show the consensus times as a function of the initial macrostate (n^jTig) for q < q c and q > q c , respectively. In the 
former case, q < q c , the expected consensus time r(nA,nB) is given by the equation f — Qt + t. A fast drop-off in 
consensus time as a function of the initial configuration, observed for q < q c [Fig. 9(a)], indicates the presence of the 
saddle point. For q > q c , there is only one stable fixed point (the absorbing one), so the consensus time starting from 
any initial configuration is short, including initial configurations in the vicinity of the previously meta-stable states 
[Fig. 9(b)]. 

Figure 10(a) and (b) present the expected number of time steps spent before absorption in each macrostate 
T(nA,ns) starting from the initial state (n q ,N — n q ) for q < q c and q > q c , respectively. From Fig. 10(a) and 
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FIG. 7: Probability of all A consensus Pa with different external influence level fa starting from macrostate (nA,riB) = 
(no, N — no) on 100- node complete graph. 

(b), according to the two peaks in each figure, the random walk before absorption spend time mainly in two areas: 
one is near the meta-stable state [close to the initial state (n g , N — n q )\, the other one is around the consensus state 
(N, 0). When q < q c , the peak around the meta-stable state is dominant in the total consensus time, while for q > q c , 
it can be ignored compared to the peak around consensus. 



We sum T(riA,nB) over these two areas separately (since the peak in Fig. 10 is very concentrated, the domain of 
summation does not matter very much), and define the time spent near consensus state as T c = J2n A >n B ^( n ^> n s) 
and that near the meta-stable state as T m = J2n A >n B T{nA,nB)- Figure 11 shows that the normalized time spent 
near consensus T c /N has the same order 0(ln(iV)) when N grows as demonstrated in the non-committed case. We 
conclude that for different q's regardless whether it is less or greater than the critical q c , the peaks around the 
consensus state have roughly the same scale and are about twice the height of the corresponding peak in the 2-word 
NG without committed agents shown in Fig. 4. Figure 12 shows that the normalized time that the random walk 
is stuck in the vicinity of the meta-stable state, T m /N, grows exponentially with N for q < q c , while it decreases 
weakly with N for q > q c . The crossover and departure from these drastically different scaling behaviors appears at 
around q = 0.08, hence our rough estimate for the critical fraction of committed agents is q c « 0.08 ±0.01. A detailed 
finite-size analysis of this crossover behavior should be performed to extract q c in the infinite system-size limit. Since 
for the consensus time we approximately have r ss T m + T c , the above findings imply that t/N ~ 0(e cN ) (where c is 
a constant) for q < q C) while t/N ~ 0(\n(N)) for q > q c . 
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FIG. 8: Probability of all B consensus 1 — Pa starting from macrostate (ua, «b) = (N/2, N/2) as a function of network size TV 
with different external influence level f's. 

VI. SUMMARY 

We studied influencing and consensus formation, in particular, the asymptotic consensus times, in stylized 
individual-based models for social dynamics on the complete graph. We accomplished this by a coarse-graining 
approach (lumping microscopic variables into macrostates) resulting in an associated random walk. We then analyzed 
first-passage times (corresponding to consensus) and times spent in each macro-state of the random- walk model. The 
method yields asymptotically exact consensus times for large but finite complete graphs of size N. Direct individual- 
based simulations can become time consuming for large systems and prohibitive even for moderately-sized systems 
when the system initially is in a meta-stable configuration, as the escape time can increase exponentially with the 
system size. The method presented here provides an alternative way to obtain the asymptotic behavior of consensus 
times, including the cases associated with extremely slow meta-stable escapes. 

After testing this framework on spontaneous opinion formation in two known models, we applied it to two scenarios 
for social influencing in a variation of the 2-word naming game. First, we considered the case when individuals are 
exposed to a global external field (or central influence). We found that the external field dominates the consensus 
in the large network-size limit. Second, we investigated the impact of committed individuals with a fixed designated 
opinion (i.e., individuals who can influence others but themselves are immune to influence). In the latter case, we 
found the existence of a tipping point, associated with the disappearance of the meta-stable state in the opinion space: 
When the fraction of committed nodes is below a critical value, consensus times increase exponentially with system 
size; on the other hand when the fraction of committed nodes is above this threshold value (tipping point) the system 
is quickly driven to consensus with weak system size dependence. 

While the method and the results presented here are applicable to the complete graph, consensus times often exhibit 
the same asymptotic scaling with the system size in large homogeneous sparse random networks [9, 11, 12], hence our 
results yield some insight how ordering and consensus can evolve in realistic social networks. In particular, one can 
better understand and predict timescales associated with reaching consensus in social networks. 
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FIG. 9: Expected normalized consensus time (r(n)/N) as a function of the initial macrostate (ua, tlb) on the complete graph 
with N = 100 nodes, (a) when the fraction of committed agents is q = 0.06 < q c ; (b) when q = 0.12 > q c . 
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FIG. 10: Expected time spent in each macro-state before consensus T(jiA,nB) on the complete graph with N = 100 nodes, 
starting from the (ua{0), ns(0)) = (n q ,N — n q ) initial macrostate, (a) for q = 0.06 < q c ; (b) for q = 0.12 > q c . 
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FIG. 11: Normalized time spent near the consensus state before consensus as a function of network size N for different fraction 
of committed agents q, including cases for both q < q c and q > q c - 
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